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ABSTRACT 

We have used the Very Large Array (VLA) at 7 mm wavelength to image five rotational transitions 
( J = 1 - 0) from three SiO isotopologues towards Orion BN/KL: ^SSiO w = 0, 1,2; and ^^SiO and 
■^"SiO V — 0. For the first time, we have mapped the ^^SiO and '^"SiO J = 1 — emission, established 
the maser nature of the emission, and confirmed association with the deeply embedded high-mass 
young stellar object commonly denoted radio Source I. The ^^SiO u = maser emission shows a 
bipolar structure that extends over ^ 700 AU along a northeast-southwest axis, and we propose that 
it traces a bipolar outflow driven by Source I. The high-brightness isotopic SiO maser emission imaged 
with a ^0''2 resolution has a more compact distribution, generally similar to that of the ^^SiO v = 1,2 
emission, and it probably traces bulk gas fiows in a region of diameter ^100 AU centered on Source I. 
On small scales of < 10 AU, however, compact ^^SiO/'^°SiO t; = and ^^SiO v — 1,2 emission features 
may be offset from one another in position and line-of-sight velocity. 

From a radiative transfer analysis based on a large velocity gradient (LVG) pumping model, we 
derive similar temperatures and densities for the optimum excitation of both ^^SiO/^°SiO v = and 
^^SiO w = 1, 2 masers, significantly higher than required for ^^SiO v — maser excitation. In order to 
account for the small-scale differences among the isotopologues {v — 0) and the main species {v = 1, 2), 
follow-up radiative transfer modeling that incorporates non-local line overlap among transitions of all 
SiO isotopic species may be required. 

Subject headings: line: profiles - masers - ISM: individual objects: Orion BN/KL - stars: formation 



1. INTRODUCTION 

The Orion Becklin-Neugeb auer/Kleinman-Low 
(BN/KL) region {D - 414 ± 7 pc; iMenten et all [20071) 
is one of only three star forming region s (SFRs) known 
to exhibit maser emissio n from ^^SiO (|Hasegawa et al.l 
119861 : IZapata et al.l 1200 8') and the only one to display 
emission from the major SiO isoto pic species comp rising 
1^0: 28SiO, ^'^SiO, ^osiO (e.g ., [GoddiiliDISool, and 
18Q. 28gji8o ( Cho et al.ll200l . 

Early single-dish and interferometric observations of 
J = 2 — 1 (A3 mm) and J = 1 — (A 7 mm) rotational 
transitions in the v = 1 and v = 2 vibrational states 
have shown similarities among spectral profiles and spa- 
tial distributions, possibly indicating that they arise in 
the same gas volumes. The vibrationally excited states 
have bee n known to be inverted since their firs t identifi- 
cation (Tfiaddiui^tin [MZI [B^hliE^ [197^ Obser- 
vations with the Very Large Array (VLA) first showed 
that the center of the v = 1 ( J = 1 — 0) SiO maser emis- 



sion coin cides with the radio continuum source denote d 
Source I (jChurchwell et al.l[l987l : iMenten fc R,eidlll995f) . 
Subarcsecond imaging at near- and mid-infrared wave- 
lengths shows no infrared counterpart to Source I, which 
sugge sts it is a deeply embedded y oung stellar object or 
YSO (|Cezarilll992l : lCreenhill et"aLl[20b4bl High-angular 
resolution imaging with Very Long Baseline Interferom- 
etry (VLBI) in both J = l- 0andJ = 2- l fines 
has shown that the v = 1 and v — 2 emission traces 
a protostellar wind and/or outflow expanding with a 
characteristic vel ocity 5, 20 km s~^ at r adii of 20-70 AU 
from Source I (iGreenhill et al.l [19981 iDoeleman et al.l 
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11999^ 'Greenhill et all l2004at IDoeleman et all 12004 
Matthews et al. 2003 Matthews et al. in prep.). 

Thermal ^*SiO v = emission is a well-known tracer 
of shocks in protostellar jets (e.g., iNisini et al.l [2007) . 
After the discovery o ft' = 0J = 2 — 1 (|Buhl et al.l 
119750 and J = 1 - (|Genzel et al.lll980[ ) fines in Orion 
BN/KL, interferometric imaging demonstrated that the 
vibrational ground state emission is part mase r and part 
thermal (J = 1-0: [Chandler fc de Pre3 [l995l: J = 2-1: 
iWright et al.][l995h . The ^SSiO v = traces a bipolar 
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structure elongated northeast-southwest (NE-SW), that 
has been interprete d either as a rotating disk ( J = 2 — 1; 
iWright et all fl99l or as a wide-angle bipolar outflow 
with characteristic sp eeds of 10-15 km s~^ {J = 1^0; 
iGreenhill e t al. 2004a\ 

The ^^SiO and ''"SiO isotopologues were first discov- 
ered towards Orion BN /KL in the v = J = 2 — 1 
line (iLovas et all 119761 : and IWolfl Il980l respectively). 
lOlofsson et alj (|198lD first proposed maser action for 
the SiO (J = 2 — 1) isotopic emission, based on ob- 
served (yearly) variability of individual Doppler com- 
ponents (^^Si O ) or t he asymmetric line shape ('^°SiO). 
iBaudry et al.l (|1998l ) confirmed maser action for the 
^^SiO J = 2 — 1 transition via high-angular resolution 
(1 — 3") imaging, but there has been no conclusive evi- 
dence of maser action for the J = 1 — transition or for 
■^°SiO in either transition. 

The physical origin of the isotopic emission is not 
yet well determined. Plateau de Bure Interferometer 
(PdBI) imaging shows that the centroids of ■^^SiO emis- 
sion (J — 2 — 1) in multiple velocity channels be- 
tween — 11 and 24 km s~^ arise ^100 A U from Source 
I, wh en observed with 1 — 3" resolution (jBaudrv et al.l 
I1998D . In contrast, the higher order J = 8 — 7 tran- 
sition of ^°SiO, observed with the Smithsonian Sub- 
millimeter Array (SMA) and comparable angular reso- 
lution, is distributed over ^ 1000 AU, elongated NE- 
SW, as i s the part-thermal, part-maser ^^SiO v = 
emission (jBeuther et al.l 120051 ) . Single-dish spectra for 
the ground-state ^^SiO and ^OSiO J = 1 - transi- 
tions also exhibit line wings that are simila r to those of 
^^SiO (±40 km s~^ ), as recently reported bv lGoddi et al.l 
()2009[ ). These are broader than the velocity extent of vi- 
brationally excited emission and invite speculation that 
they might trace moderately high-velocity components 
of the flows around Source I fcToddi ct al. 2009) . 

Sub-arcsec resolution imaging is necessary to establish 
unambiguously the location and the nature (maser vs 
thermal) of the isotopic emission within the BN/KL re- 
gion, in particular where the high- velocity isotopic emis- 
sion arises within the SiO flow. With this purpose, we 
conducted simultaneous observations of the v = 0,v = 1, 
and i; = 2J=l — transitions of ^*SiO and the 
w = 0J = l — line of the ^^SiO and '^'^SiO isotopologues 
at 7 mm wavelength with the VLA. For the first time, we 
have mapped the ^^SiO and the ■^"SiO w = 0J = l — 
lines, measured the proximity of the emission regions to 
Source I, and established the presence of population in- 
version. Simultaneous observation of SiO masers from 
three vibrationally excited transitions and three isotopo- 
logues, provide constraints on SiO excitation conditions 
in the surrounding medium. 

The observational setup and data calibration proce- 
dures are described in §2. Spectral profiles and images of 
different maser transitions are presented in §3. In §4, we 
discuss the gas kinematics probed by each species (§4.1), 
the presence of maser action from ^^SiO and ^"SiO (§4.2), 
and the excitation conditions derived from preliminary 
radiative transfer modeling (§4.3). Finally, conclusions 
are drawn in §5. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 



TABLE 1 
Parameters of Observations. 



Program 


Date 


VLA(i) 


Transitions 


^rcst 


T{2) 






Conf. 


J = 1 - 


- 


(MHz) 


(hr) 


AG776 


2008/02/02 


B {13V) 


28 SiO V 


= 


43423.858 


2.3 


AG776 


2008/02/02 


B {13V) 


28 SiO V 


= 1 


43122.080 


4.6 


AG776 


2008/02/02 


B {13E) 


28 SiO V 


= 2 


42820.587 


4.6 


AG776 


2008/02/02 


B {13V) 


29 SiO V 


= 


42879.846 


2.3 


AG776 


2008/02/02 


B {13E) 


30 SiO V 


= 


42373.340 


4.6 


AG575 


2000/01/01 


B {19} 


28 SiO V 


= 1 


43122.080 


0.5 


AG575 


2000/01/01 


B {19} 


29 SiO V 


= 


42879.846 


0.5 


AG575 


1999/08/28 


A (10) 


28 SiO V 


= 1 


43122.080 


0.8 


AG575 


1999/08/28 


A (10) 


29 SiO V 


= 


42879.846 


0.4 



Note. — 

^ Spacing of antennas: B-configuration has baselines of 0.1 to 11 km. A- 
configuration has basehncs of 0.5 to 35 km. The number of antennas partici- 
pating is in parentheses, trailed by "V" to indicate VLA and "E" to indicate 
EVLA. Prior to 2008, the number of participating antennas corresponds to 
the number equipped with high-frequency receivers. All setups placed ei- 
ther ?; — 1 or f — 2 ^^SiO in one IF to enable high-accuracy calibration. 
^ Approximate on-source integration time. 

Observations of the SiO emission from Source I were 
conducted at 7 mm using the Very Large Array (VLA) 
of the National Radio Astronomy Observatory (NRAO)^. 
We present data obtained on 2008 February 2, along with 
archival data from 1999 August 28 and 2000 January 1. 
The data from 2000 and 2008 were obtained while the 
VLA was in the B-configuration, yielding a resolution of 
~0"2. The data from 1999 arc from the A-configuration, 
giving a resolution of ~0''5. Table [T] summarizes the 
observational parameters of the different programs. The 
following description refers to the program observed in 
2008. 

We observed a total of five J = 1 — rotational tran- 
sitions from three different SiO isotopologues: ^^SiO 
V = 0,1,2; 29siO V = 0; and ^OSiO v = (Table [T]). 
To achieve high-accuracy astrometry, we cycled our ob- 
servations among the various transitions, observing two 
at a time. We included emission from the w = 1 or u = 2 
states of ^*SiO in each of the setups because these were 
strong enough to be used to calibrate the troposphere 
and system gain, thus facilitating detection of the weaker 
emission from the other isotopologues. At the time of 
the 2008 observations the VLA included 13 EVLA^ an- 
tennas that were able to simultaneously observe two sep- 
arate bandpasses (IFs) up to 2 GHz apart. In contrast, 
the maximum separation between IFs allowable with the 
remaining (VLA) antennas was ~ 400 MHz, depending 
on details of the available local oscillator settings. In 
order to optimize the use of the VLA and EVLA an- 
tennas, respectively, we therefore divided the array into 
two sub-arrays of 13 antennas each. The VLA antennas 
observed the ^^SiO ?; = 1 maser at 43122 MHz m one 
IF (single polarization) and the ^^SiO and ^^SiO v = 
fines at 43424/42880 MHz alternately in the other IF 
(single polarization). The EVLA antennas observed the 
^^SiO w = 2 maser at 42821 MHz simultaneously with the 

^ NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

2 The Expanded Very Large Array is a joint USA, Mexico, and 
Canada project to greatly improve the observing capabilities of 
the VLA array by installing new electronics {fiber optic cables, 
large bandwidth receivers, new correlator, etc.) while keeping the 
existing infrastructure (antennas, railroad tracks, etc.). 
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^°SiO t; = line at 42373 MHz (all single polarization). 
For 30SiO w = 0, a tuning offset of 2.83 MHz was intro- 
duced alternately with the main tuning to enable a search 
for possible high- velocity emission (up to 70 km s~^), 
who se presence was suggested in earlier single-dish spec- 
tra (|Goddi et al. 2009 ). No such emission was detected 
by the present observations. The total on-source inte- 
gration times for the various SiO transitions are given in 
Table [l] and the frequency pairings used in the observing 
programs are summarized in Tabled) The ^®SiO, ^^SiO, 
and ^''SiO v = transitions were each observed with 
a 12.5 MHz total bandwidth (88kms~^) with a spectral 
resolution of 390.6 kHz (2.7kms-i). The ^^SiO w = 1,2 
emission was observed with a 6.25 MHz bandwidth and 
97.656kHz (0.68 kms-^) channels (Table[2l). 

We added several continuum-mode setups to the 2008 
observations to estimate any frequency dependent (i.e., 
delay) terms to the astrometry error budget that might 
affect the relative astrometry between the pairings of 
transitions observed in spectral-line mode. In particu- 
lar, we performed observations of the respective ^^SiO 
maser {v — 1 for the VLA subarray, v = 2 for the EVLA 
subarray) in one IF (single channel, 6.25 MHz band- 
width, dual polarization) together with observations in 
a 50 MHz continuum channel in the other IF (also dual 
polarization), offset in frequency by ±350 MHz from the 
maser. The goal was to detect and image the BN object 
and use its position to align the various SiO emissions. 
Unfortunately these observations on the EVLA subarray 
failed, so we have aligned the v — 1 and v = 2 emis- 
sions with respect to each other using other techniques 
(Section [2Jl). 

We derived absolute flux calibration from observations 
of 3C 286 {Fi, — 1.5 Jy), and bandpass cahbration from 
observations of 3C 84. The phase offsets between the 
paired IFs were tracked through observations of QSO 
J0541-0541 every 25 minutes. The strong ^^SiO v = 1 
and V = 2 masers (3> 100 Jy) were used to track tro- 
pospheric phase fluctuations on 10-s timescales, as de- 
scribed in Section [2T2l All spectral-line observations were 
observed in a fixed-frequency mode, and Doppler correc- 
tions were applied during the data reduction. 

Observations of ^^SiO v = and ^^SiO w = 1 in 
1999 and 2000 (Table [1]) were made with setups simi- 
lar to those used in 2008 for the VLA antennas. We 
adopted a 2-IF, single polarization mode, with band- 
widths of 6.25 MHz. Monitoring of phase offsets between 
the IFs was accomplished with periodic observations of 
QSO J0541— 0541, and bandpass calibration relied upon 
observations of QSO J0530+135. 

2.2. Calibration procedures 

Calibration and data reduction were carried out us- 
ing the Astronomical Image Processing System (AIPS). 
First, the phase offset between the two simultaneously- 
observed IFs due to the electronics was removed, so that 
subsequent tropospheric calibration from the ^'^SiO v — 1 
and V — 2 transitions could be apphed directly to the IF 
containing the weaker line. The phase offset was deter- 
mined from observations of J0541— 0541 in the follow- 
ing manner. Using the 12.5 MHz bandwidth IF, calibra- 
tion solutions were determined for the short-term tropo- 
spheric phase fluctuations on J0541— 0541 and applied to 
J0541-0541 in the 6.25 MHz IF. A subsequent phase cal- 



ibration on the 6.25 MHz IF using J0541-0541 and av- 
eraging over scan lengths (^ 2 min) then gave the phase 
offset of this IF relative to the 12.5 MHz, and was ap- 
plied to all sources in the dataset. Amplitude calibra- 
tion was subsequently derived from J0541— 0541, after 
correcting for the short-term tropospheric phase fluctu- 
ations. We also derived separately phase solutions from 
2-minute scan averages of J0541— 0541 and applied them 
to Source I. Corrections for the shape of the antenna- 
based bandpasses were also applied, along with Doppler 
corrections. 

Finally, corrections for fluctuations in the tropospheric 
and instrumental phase and amplitude were derived ev- 
ery 10 seconds from channel averages of the strong ^^SiO 
masers (Vlsr = —7.2 to —3.1 km s^^ for the v = 1 line, 
and —9.4 to —5.3 km s"^ for the v = 2 line) using self- 
calibration and applied to the IF containing the weaker 
line. As a consequence, maps of the ground-state tran- 
sitions are (independently) referenced to the positions of 
the ^*SiO w = 1 or ti = 2 emission. 

Images were made of Source I using the AIPS task 
IMAGR with a "robust" parameter of 0, resulting in a 
synthesized beam size of '^0"2 for the B configuration 
data and ~0"05 for the A configuration data. The images 
were 1024 pixels on a side with cell sizes 0"02 (2008) 
or 0"005 (1999) and covered fields-of-view of 20" and 
5", respectively. For the ^^SiO images the sensitivity is 
dynamic range limited by the sidelobes of the strongest 
emission peaks and thus varies among different frequency 
channels and different transitions. Table [3] summarizes 
the imaging results. 

2.3. Positional registration among different transitions 

Phase and amplitude self-calibration enable coherent 
integration and imaging but residual astrometric errors 
remain, due to the ^ 1° separation between calibrator 
(J0541-0541) and target. To remove frequency and sky 
position dependent errors (i.e., the residual delay error 
due to the atmosphere), we used as a position refer- 
ence for registration the peak of the continuum emis- 
sion from BN, which is the strongest in Orion BN/KL 
{F,y ^ 25 mjy at 7 mm) and only ~ 10" from Source 
I. For that purpose, we calibrated the continuum-mode 
dataset from 2008 following a procedure similar to that 
described in Sect. 12.21 for the line-mode dataset. The 
band-averaged signal-to-noise ratio (SNR) for the v = 1 
and V = 2 maser emission over 6.25 MHz was high enough 
that the narrow band data enabled coherent integration 
and detection of BN in the broad band. 

By directly comparing the positions of BN in bands 
tuned to the v — 1 emission and offset by ±350 MHz, 
we established the angular offset per MHz introduced 
by calibrating one band by another: Aa = 0.007 ± 
0.003 mas/MHz and AS = -0.016 ± 0.004 mas/MHz. 
Unfortunately, the EVLA subarray tuned to the v = 2 
emission did not provide usable data in the continuum 
mode setups. We used the offset derived with the VLA 
subarray to correct the positions of all pairs of transitions 
with 1 mas accuracy. 

Since the v = 1 and v = 2 transitions served as 
calibration references for different subarrays and were 
self-calibrated independently, without an EVLA detec- 
tion of BN it is necessary to register them separately to 
each other. For that purpose we measured the intrinsic 
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TABLE 2 

Frequency pairs in the observing setup 



Prog. 


Mode 


IFA(-C) 


BW 


CW 


IFB(-D) 




BW 


CW 






(MHz) 


(kHz) 




(MHz) 


(kHz) 


AG776 


Line 


28SiO v = 


12.5 


390.6 


28SiO V = 


1 


6.25 


97.7 


AG776 


Line 


29siO v = 


12.5 


390.6 


28SiO V = 


1 


6.25 


97.7 


AG776 


Line 


sosio „ = 


12.5 


390.6 


28SiO V = 


2 


6.25 


97.7 


AG776 


Cont. 


= 1 + 350 


50.0 




28SiO V = 


1 


6.25 




AG776 


Cont. 


I, = 1 - 350 


50.0 




28SiO V = 


1 


6.25 




AG776 


Cont. 


V = 2 + 350 


50.0 




28SiO V = 


2 


6.25 




AG776 


Cont. 


■u = 2 - 350 


50.0 




28SiO V = 


2 


6.25 




AG575 


Line 


29SiO v = 


6.25 


97.7 


28SiO V = 


1 


6.25 


97.7 



Note. — Cont. — Continuum and Line — Spcctral-linc mode: 64 (32) channels of 97.7 
(390.6) KHz width (CW) within a 6.25 (12.5) MHz bandwidth (BW). 



TABLE 3 
Summary of Observational Results 



Date 


Transitions 
J = 1 - 




Synthesized Beam'^' 

eAf(") xe™("); P-A.n 


RMS(2) 
(mjy beam~^) 


Blue Peak(^) 
V(km s-i) I(Jy beam-l) 


Red Peak(^) 
V(kms-i) I(Jybeam-i) 


AG776 


28SiO V = 





0.23 X 0.19; 34 


7-113 


-5.8 


10 


13.1 


4.9 




28SiO V = 


1 


0.23 X 0.19; 33 


15-170 


-4.5 


1590 


14.5 


930 




28SiO V = 


2 


0.24 X 0.15; 20 


6-40 


-6.6 


470 


20.0 


380 




29SiO V = 





0.22 X 0.19; 25 


5 


-3.7 


0.044 


15.5 


0.12 




30SiO V = 





0.25 X 0.15; 21 


5 


-4.6 


0.068 


14.8 


0.072 


AG575 


28SiO V = 


1 


0.20 X 0.17; 35 


11-2003 


-3.2 


340 


14.5 


210 




29SiO V = 





0.20 X 0.17; 7 


7-27 




< 0.05 


16.2 


0.34 


AG575 


28SiO V = 


1 


0.061 X 0.045; 53 


50-200 


-3.9 


600 


13.8 


400 




29SiO V = 





0.059 X 0.042; 45 


16 




< 0.07 


15.9 


0.840 



Note. — 

' The synthesized beam results from a robust weighting. 

' 1(7 RMS noise in a 0.7 or 2.7 km channel, based on histogram fits to the pixels over the channel image. 
' Velocity and intensity of the emission peak on the blue and red side of the spectra. 



structure of the v — 1,2 emission using VLB A images 
from a monthly monito ring of Source I over several years 
(jMatthews et al.l[2007t Matthews et al. in prep.). After 
averaging in velocity to account for different spectral res- 
olution, the VLBA images were convolved with the VLA 
beam, and the offset between the v — 1 and v — 2 emis- 
sion was measured. We did tests using 4 to 5 different 
velocity ranges in the v = 1,2 spectra at 4 to 5 dif- 
ferent epochs from the VLBA monitoring and obtained 
consistent results with an offset of Aa = 14 ± 3 and 
A5 = 0.5 ± 3 mas. The accuracy in the computed offsets 
is dominated by the variation of the structure and of the 
relative locations of the v = 1 and v = 2 maser centroids 
from channel to channel and from epoch to epoch, rather 
than by the uncertainties in measuring the locations of 
the V ^ 1,2 emission peaks. The quoted uncertainty is 
then the variance of the offset values measured in several 
channels and at different epochs. 

After applying the angular offset per MHz relation de- 
rived above and the offset between the v = 1 and v = 2 
emission, we estimate that the maps of the various tran- 
sitions are aligned on the sky with accuracy of 3 mas in 
each direction. 

2.4. Absolute astrometry 

The application of the phase solutions derived from 2- 
minute scan averages of J0541— 0541 to the maser data 
imposes an absolute astrometry (Section 12. 2p . How- 
ever, because J0541— 0541 was observed relatively infre- 



quently (every 25 minutes) compared with the timescale 
of tropospheric and instrumental phase variations, the 
absolute astrometry is poor, and further alignment is 
necessary. In this case, absolute positions were es- 
tablished by comparing posi tions of BN i n our maps 
with the position derived by iGomez et al.l ()2008f ) from 
VLA observations conducted in May 2006 at 8.4 GHz, 
a(J2000) = 05''35'"14?1099, 5{J2QQQ) = -05°22'22"741 
(positional accuracy ~0"01) and their estimated proper 
motion (14 = —5.3 mas yr~^, Vy = 9.4 mas yr""'^). For 
2008 February, we obtain a(J2000) = 05'*35"14f 1093, 
(5( J2000) = -05°22'22"7253. The offset between the po- 
sition of BN measured in our maps and and t he pr oper- 
motion-corrected position from iGomez et al.l ()2008f ) was 
applied to the SiO images, whose absolute position is 
then measured to an accuracy of 0''01. 

2.5. Fitting positional centroids of maser emission 

The emission in each velocity channel is unresolved in 
all transitions except the ^^SiO u = 0. For the latter, in 
order to establish the morphology and kinematics of the 
extended emission, we produced an intensity map inte- 
grated over the emission velocity range. For each unre- 
solved emission component, a two-dimensional ellipsoidal 
Gaussian model is fitted in a 100 x 100 pixels area of each 
channel map. Since the angular extent of the maser emis- 
sion in each velocity channel is smaller than the synthe- 
sized beams of the VLA in the observing configurations, 
spatially separated maser components at similar veloci- 
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ties are blended together. Hence, the Gaussian-fitted po- 
sitions actually locate the intensity-weighted centroids of 
SiO emission in each frequency channel. In the follow- 
ing, we will refer to these emission centroids in individual 
velocity channels as maser features. 

The relative positional errors of maser features are lim- 
ited by the SNR of the channel maps and by systematic 
errors in the calibration of the bandpass. In the case 
of weak signals (e.g., ^^SiO and ■^"SiO), the positional 
uncertainties are noise limited, and are given by SO = 
0.56'b/SNR, where 63 is the FWHM of the synthesized 
beam, and SNR is the peak intensity divided by the RMS 
noise in a particular velocity channel. For the 2008 maps 
{9b =0''2), the relative positional uncertainties are typi- 
cally ~ 10 mas for ^^SiO and ^"SiO (-50 mJy beam"! 
peak and ^^5 mJy beam^^ RMS), whereas for the 
1999 maps {9b 50 mas), the typical uncertainties 
are ^ 1 mas for -^^SiO (~500 mJy beam~^ peak and 
^17 mJy beam^^ RMS). For the v — 1,2 emission, 
the positional uncertainties are dominated by frequency- 
dependent errors in calibrating the bandpass. From the 
continuum spectrum of 3C 84, we measure a channel to 
channel phase noise < 1°, which corresponds to an un- 
certainty of < 0.6 mas for a 0"2 synthesized beam. 

3. RESULTS 

We imaged the ^^SiO v = 0,1,2 and the ^^SiO and 
30siO w = 0J=l-0 transitions observed in 2008 (0"2 
beam), as well as the ^^SiO v = 1 and ^^SiO v — 
transitions observed in 2000 (0''2 beam) and 1999 (0"05 
beam). 

In each epoch and for each transition, we produced 
spectra by mapping each spectral channel (0.68 or 
2.7 km s~^ wide) and then summing up the flux density 
in each channel map. Figure [1] compares for each tran- 
sition the integrated flux densities of the VLA channel 
maps at different epochs with the G ET single-dish spec - 
tra (observed in 2007 November 10: iGoddi et al]l2009f ). 
Multiple transitions show similar double-peaked profiles 
and velocity extents: blueshifted emission in the velocity 
range from —10 to km , redshifted emission in the 
interval from 10 to 25 km s""'^ , and systemic emission 
in the range from to 10 km s^^ . Comparison of the 
GET (2007 November) and the VLA (2008 February) 
spectra, shows an increase in the flux density of about 
100% in the emission oi v = 1 (both blue and red com- 
ponents) and V = 2 (only red component), whereas the 
^^SiO V — emission does not show any significant vari- 
ation. Conversely, less than — 30% of the single-dish 
emission has been detected in the VLA channel maps 
for ^^SiO and ■^"SiO. The single-dish spectral profiles of 
the latter two transitions appear much broader than the 
VLA profiles and show high-velocity wings undetected 
in the VLA channel maps (Fig. [T]). Interestingly, the 
^^SiO spectrum has undergone significant changes over 
8-9 years (the separation between VLA epochs). In the 
spectra taken in 1999 and 2000, the redshifted compo- 
nent is 10 and 5 times stronger, respectively, than in the 
2008 spectrum while no emission from the blue compo- 
nent is detected above a 4tT ~ 50 — 70 mJy noise level 
(Fig.©. 

The velocity integrated intensity map of ^^SiO v = 
emission consists of two bright peaks offset — ±0"5 
from Source I (Fig. [2l left panel). The bipolar structure 




"SiO v=l Z 

VIA (2008/02/02)- 

VIA (2000/01/01): 

VIA (1999/08/28): 

GET (2007/11/09)- 



60 

""siO v=2 - 

VLA (2008/02/02)^ 
GBT (2007/11/09)] 

^"siO v=0 ^ 

VIA (2008/02/02): 
GBT (2007/11/09): 

'"'siO v=0 ; 

VLA (2008/02/02): 
VLA (2000/01/01): 
VLA (1999/08/28)^ 
GBT (2007/11/09): 



VLA (2008/02/02)- 
GBT (2007/11/09): 



CH3OCH3 



20 

Velocity (kin/s) 



Fig. 1. — Spectra of various SiO transitions observed toward 
Orio n BN/KL with the VLA during several epochs. GBT spectra 
from IGoddi et aLI 12009) are also shown for comparison. In 2008, 
the VLA spectral resolution was 97.7 kHz (0.68 km ) for the 
28SiO D = l,2J=l-0 lines and 390.6 KHz (2.7 km s"! ) for 
the 28siO, 29siO, and ^OSiO i; = OJ = l- lines. In 1999 
and 2000, SiO v = 1 and ^'^SiO v = were observed with a 
0.68 km s~^ resolution. The large differences between GBT and 
VLA V = 1,2 spectra is likely to be due to the time separation 
of observations and the strongly maser-nature and compactness 
of the emission. In the case of the isotopologues, most probably 
the observed difference is due to an extended, thermal component 
detected by the single-dish but resolved out by the interferometer. 
Note that the same flux density scale is adopted for ^^SiO and 
30SiO. 

extends over ~ 700 AU along a NE-SW direction and 
is si milar to previous B IMA images of the J = 2 — 1 
fine (Wrig ht et allll995f) . Both lobes emit in a similar 
velocity interval, from -9 to 16 km s^"'^ . 

The emission centroids from ^*SiO v ~ 1,2 and 
^^SiO and ^°SiO v = are clearly located at the cen- 
ter of the 2* SiO w = bipolar structure (Fig. [H left 
panel) and are distributed in a region of size 0"2 around 
Source I (Fig. [21 right panel). They show generally simi- 
lar spatial and velocity distributions and can be divided 
into three groups: blue features detected southeast (SE) 
of Source I include emission with line-of-sight (L.O.S.) 
velocities blueshifted with respect to the systemic LSR 
velocity of the region (5 km s""'^), red features located 
to the north-west (NW) have emission with redshifted 
L.O.S. velocities, and systemic features, which are weaker 
and located in between the two other groups, are com- 
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paratively close to the systemic velocity. However, some 
differences can be noted among various transitions. Com- 
pared tov = 1, the ^*SiO v = 2 emission appears located 
in a slightly inner layer, closer to Source I, even if there is 
a significant overlap. In addition, the systemic features 
appear stronger (relative to blue and red peaks) for v = 2 
than for v = I. This is consistent with the findings of 
higher-resolution imaging of the ^^SiO v — 1,2 emission 
with the VLBA (jMatthews et al.ll2007t Matthews et al. 
in prep.). 

As far as the isotopologues are concerned, although 
correspondence among ground-state isotopologue and vi- 
brationally excited emission distributions is generally 
good, a one-to-one correspondence in positions and 
L.O.S. velocities in features among different transitions 
is clearly absent. In addition, the brightest features from 
the main species and isotopologues are not co-spatial 
(with the exception of the ^^SiO red features in the 
NW area). Finally, while most of the blue and systemic 
^^SiO features also have corresponding '^°SiO features, 
this is not necessarily true for the red features. 

We caution that, owing to low angular and spectral 
resolution, blending may hamper detailed comparisons 
of emission centroid positions. For instance, higher an- 
gular (0''05 beam) and spectral (0.68 km s^^ ) resolu- 
tion imaging with the VLA reveals that the ^^SiO v — 1 
maser emission does not arise from two (blue and red) 
arcs (as shown in Fig. but resolves into an "X"- 
shaped pattern, consisting of four "arms" , extending to- 
wards north (N), west (W), south (S), and east (E), 
plus a bridge of emission joining the S and W arms of 
the "X" (Figure [3]) . The X-shaped structure is consis- 
tent with sub-milliarcsecond resolution VLB A imaging 
(jGreenhill et al.lll998t iDoeleman et "all 12004 Matthews 
et al. in prep.). In the high-resolution VLA map, most 
of the red features arise from equally bright components 
in the N and W arms, the blue features from components 
in the S and E arms, and the systemic features from the 
bridge (Figure [3]). As a consequence, with a 0"2 angular 
resolution and 2.7 km s"'^ spectral resolution, centroids 
of both red and blue features fall between N and W arms 
and S and E arms, respectively, resulting in the two (red 
and blue) arcs seen in Fig. [2l Similar arc morphologies 
have been previously seen in low -resolution imaging at 
7 mm with the VLA (0"2 beam - |Menten & Reid''1995^ 
and at 3 mm with BI MA (0!^5 beam - Wright et al. 1995) 
and PdBI (2" beam - iBaudrv et all fl998). 

The high resolution VLA image allows accurate local- 
ization of the ^^SiO emission towards the N arm of the 
X traced by the v = 1 emission (Fig. [3]), which invites 
speculation that the presence of isotopologue emission 
within the arcs on either side of Source I (Fig. when 
observed at lower resolution, is indicative of an origin in 
the four arms. Nonetheless, detailed comparison requires 
caution in light of the apparent 2 — 10 AU offset between 
^^SiO and ^^SiO features of similar velocity (Fig. [3]). 

4. DISCUSSION 

4.1. Gas kinematics around Source I 

The NE/SW-oriented bipolar structure traced by the 
^^SiO vibrational ground state has been previously inter- 
preted either as a rotating disk, base d on maps of the J = 
2 — 1 emission at 3 mm wavelength ()Wright et al.lll995D . 
or as a bipolar outflow, based on the distributions of 



hotspo ts of J = 1 — emission at 7 mm (jGreenhill et al.l 
I2004a^ a nd the distribution o f the J = 8 - 7 emission at 
0.9 mm (jBeuther et al.l [20051) . The data presented here 
are consistent with the outflow hypothesis, where we cite: 
(i) a fully develope d bipolar geometry whe n imaging de- 
tects both hotspots (jGreenhill et al.ll2004af ) and extended 
structure (this work); (ii) lobes that exhibit similar mixes 
of blue and redshifted components, indicative of outflow 
close to the plane of the sky rather than rotation; and 
(iii) the putative outflow in the vicinity of Source I is 
aligned with the principle axis of the "18 km s~^ " out- 
flow in the BN/KL region (jGenzel et al.lll981[ ). which is 
suggestive of Source I as a possible contributor. 

The X-shaped distribution of vibrationally excited 
^^SiO emission has been interpreted as tracing ci ther a 
bicon i cal outflow with a SE-NW axis (Grccnhi ll et al.l 
Il998t IDoeleman eTall [1991), funnel-like flows arising 
from the surface of a disk with a NE-SW rotation 
axis (jGreenhill et al.l [2004af) . or the interface between 
infaU ing and outflowing material ([Cunningham et al.l 
[2005h . 

Emission from both vibrationally excited states of the 
main species and the ground state of the isotopologues 
traces the innermost (;$100 AU) circumstellar gas around 
Source I (Fig. [2j and [3]). Notwithstanding similarities, 
it is not however possible to make an exact position- 
velocity- intensity pairing of the ^^SiO v — 1,2 and 
29SiO/30SiO w = features (Sect. [3|). Hence, it re- 
mains possible that the main and isotopic species might 
be excited in somewhat different volumes and in prin- 
ciple might require slightly different physical conditions 
(see Sections 14. 3p . A correspo nding uncertainty in in- 
terpretation has been raised by iBaudrv et al.l ([19981 ) to 
explain differences in position between the ^^SiO v = 1 
and ^^SiO w = 0J = 2 — 1 emission in 3 mm wavelength 
PdBI images. 

Previous imaging of thermal emission from ^^SiO and 
^°SiO (u = J = 8 — 7) in the submm-wave band with 
the SMA has shown that both ^^SiO and ^"SiO emis- 
sions subtend similar angles on the sky (size~ 1000 AU), 
are elong ated NE-SW, a.nd co ver about the same veloc- 
ity range ([Beuther et al.|[2b05f ). The natural inference is 
that both species may trace inner reaches of the NE-SW 
outflow. Our VLA maps, however, do not show any in- 
dication of isotopic emission within the NE-SW outflow 
traced by ^*SiO v = emission. The actual spatial struc- 
ture of the ^^SiO and ■^"SiO u = OJ=l — emission 
might be more complex and extended than illustrated in 
Figures [2j and [3| The large difference in flux between the 
single-dish and the interferometer spectra (Fig. [T]), may 
be in fact ascribed to an extended, low-surface bright- 
ness component, detected by the GBT but filtered out 
by the VLA. Indeed, in our VLA maps we arc insensitive 
to emission with an integrated flux density < 100 mJy 
and that extends over > 1" (note that the line wings 
observed in the GBT spectra for ^^SiO and '^"SiO are 
'--^100 mJy). Hence, the ^^SiO and '^"SiO emission im- 
aged with the VLA traces only the highest brightness re- 
gions in the proximity of Source I. With the present VLA 
data, we cannot establish the location of the moderately 
high- velocity component of the protostellar outflow from 
Source I, as traced by the v = wing emission observed 
with the GBT. 
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Fig. 2. — Distribution of isotopic SiO maser emission in Orion BN/KL observed on 2008 February 2. (Left) Emission map of ^**SiO v = 
J = 1 — is integrated from -11 to 21 km s~^ . The contour levels are 1, 5, 10, 20, 30, 40, 50, 60 Jy b eam~^ km s~^ . The map shows a 
bow tie shape similar to that observed previously at 3 mm wavelength with BIMA l lWright et aLllTQQSI) . Source I is located at the center 
of the ^^SiO V = flow, where the emission from other isotopologucs and vibrational states also lie. The beam (0"22 X 0"19, P.A. 25°) is 
shown on the right upper corner. (Right) Emission centroids of ^*SiO d = 1 (filled circles) and v = 2 (crosses), ^^SiO (filled triangles) and 
•^''SiO D = (open squares) as a function of velocity superimposed on the velocity integrated intensity map of ^*SiO v = (contours). Color 
denotes L.O.S. velocity (color scale on the right-hand side). The sizes of circles and crosses scale logarithmically with the peak intensity 
of the V = 1,2 emission. For ■^^SiO and '^"SiO the position uncertainties reflect formal errors in model flts. The adopted absolute position 
of So urce I is a(J2000) = 05''35'"14f 5146 ± Of 0007, <5(J2000) = -05°22'30!'5655±0;'01, after correction for proper motions l IGomez et al.l 
[2003 ): the size of the star indicates uncertainty. The linear and angular spatial scales are given in the lower left and upper right corners 
of the panel, respectively. The relative alignment between pairs ^*SiO v = 1 and ■^*SiO/^^SiO d = or ^*SiO v = 2 and ^"^SiO j; = is 
accurate to ~3 mas. The uncertainty in the absolute position is ^^10 mas. 



4.2. Nature of the isotopologue emission 

The 29 SiO and 30 SiO w = J = 1 - emissions 
have been imaged at sub-arcsecond resolution for the 
first time. [Chandler fc de Pred (ligPS *) showed that the 
2^SiO ?; = J = 1 - emission in Ori on BN/KL is 
art t hermal and part maser {Tb - 10^ K). lBaudrv et all 
1998f ) reached the same conclusion for the ^^SiO v = 
J = 2 - 1 emission {Tb - 10^ K). 

On one hand, the broad profiles and the large flux den- 
sities at systemic velocities in the ground-state emission 
as compared to the vibrationally excited states (Fig. [1]) 
clearly indicates thermal emission from all species. On 
the other hand, the high flux density measured for the 
strongest features in the interferometric maps suggests, 
similarly to ^^SiO, non-thermal emission for ^^SiO and 
3° SiO w = OJ=l — Oas well. Although our ob- 
servations do not resolve the emission in individual ve- 
locity channels (and hence we do not know the actual 
size of individual features), we may use the synthesized 
beamwidth and the observed peak flux density at differ- 
ent angular resolutions to estimate a minimum brightness 
temperature. For ^^SiO, the peak intensity measured 
in the 16 km channel map with a 0''2 resolution is 
0.12 Jy beam~^, which corresponds to a line brightness 
temperature of > 1900 K. In the 1999 observations, with 
an angular resolution of 0"05, the same velocity compo- 
nent has a peak intensity of ^ 0.84 Jy/beam, correspond- 
ing to a brightness temperature of > 2 x 10^ K. Hence, 
the bright feature at 16 km s~^ is undoubtedly a maser. 



Temporal variations in the ^^SiO line profile (Fig. [T]) fur- 
ther support non-thermal processes in the excitation of 

29SiO. 

The '^'^SiO emission is likely maser in nature as well. 
For the strongest '^"SiO emission in the -4.8, -1.8, and 
14.8 km s^^ channels, we infer a brightness temperature 
of > 1200 K (for a 0"2 beam), which is comparable to 
the value inferred for ^^SiO from maps with the same 
resolution. The estimated value is well above the ki- 
netic temperature observed even in the ne arby hot core 
(150-300 K) fe.g.. ' Cxenzel fc Stutzkil [l98?) . but it could 
be plausible in principle in shocked gas, where the SiO 
molecule is believed to be produced and excited. In par- 
ticular, the SiO molecule is produced in C-type shocks by 
the injection into the gas-phase of Si by grain sputtering 
and/or grain-grain collisions, followed by gas- phase reac- 
tions with O/O2 (.Schilke et a l. 1997; C aselh et al. 1997). 
At typical speeds of C-type shocks (10-50 km s~^ ), the 
temperature of the post-shocked gas rises up to 1000 K 
and then progressively decreases to its pre-shock values. 
Multi-line SiO observations in protostellar jets give ki- 
netic temp eratures < 1000 K (typical values in t he range 
100-500 K: [Nisini et al.|[2007HCal3rit et al.ll2007D . Hence, 
if the emission were thermal, the inferred brightness tem- 
perature would be only marginally consistent with the 
temperatures that could be achieved via C-type shocks 
and/or with typically observed kinetic temperatures in 
SiO protostellar jets. Based on similarities (e.g., com- 
pactness) with all other masers lines (^^SiO v = 1,2 and 
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Fig. 3. — Centroids of ^*SiO v = I [filled circles) and 
■^^SiO V = {filled triangles with error bars) compact emission ob- 
served on 1999 August 28. The ^SSiO v = l emission distribution 
traces a 4-arm structure, redshifted in the north and west arms and 
blueshifted in the south and east arms. Note also a weaker bridge 
of emission connecting west and south arms. The ^^SiO emission 
is concentrated towards the north arm {dashed rectangular area). 
The ellipse on the top left indicates the synthesized beam (0"06 
X 0"04; P.A. 53°). Since the beam size is comparable to the arm 
extensions, the v = 1 emission from each arm is not resolved in the 
map, and the specific structure within each arm may be artifactual 
owing to spatial blen ding (compare with the VLBA images from 
IGreenhill etaI|[T99^ . Note that the same physical scale is used in 
both Fig. |2] and Fig.|3] 

^^SiO V = 0), the "^^SiO emission is more likely to be a 
maser. 

For some of the weak ^^SiO and '^"SiO features we ob- 
tain brightness temperatures in the range 260-400 K, in- 
dicating that the isotopic v = emission imaged with the 
VLA is probably part thermal and part maser as for the 
^^SiO V — emission. The difference between the single- 
dish and interferometric spectral profiles indicates that 
the VLA in B-configuration does not have the brightness 
temperature sensitivity (RMS noise Tj, ~80 K) to detect 
the weaker and extended thermal component fSect. |4?T|) . 

Based on our single-dish spectra and interferometric 
images, we conclude that the v = emission from all 
species is a combination of thermal and maser compo- 
nents. 

4.3. Physical conditions from a radiative transfer 
modeling 

The ground-state isotopologue and the vibrationally 
excited state emissions appear to arise from gas within 
~ 100 AU of Source I, but possibly at somewhat different 
locations within the volume. We have undertaken pre- 
liminary radiative transfer calculations to assess whether 
putative offsets might be due to excitation requirements 
for the different transitions. 



We used a radiative transfer code based on th e Larg e 
Velocity Gradient (LVG) model of iDoel et all (|1995^ . 
that includes both coUisional and radiative pumping but 
does not treat effects of maser saturation or line over- 
lap. We carried out model calculations using 600 en- 
ergy levels = — 4; J = — 39 for each isotopic 
species) . The main input parameters for the LVG model 
are reported in Table [H The model background con- 
tinuum assumes a blackbody function of = 10^ K. 
The radiation field is reduced by a geometric dilution 
factor, W — l/4(i?^3^(j/i?^as), where i?i.ad is the radius 
of the radiative source and i?mas is the distance to the 
maser region. The lower limit of i?rad corresponds to 
a B star photospheric radius (7 x Rq) (jZombeckll2007f) 
and the uppe r limit roughly ap proximates the radio con- 
tinuum size (jReid et al.ll2007f ). We do not take into 
account emission from dust around Source I as dust is 
likely to be sublimated at radii < 100 AU, where high- 
density SiO in gas phase is present. A velocity gradient 
dV/dr = 0.5 km s~^ AU~^ is assumed based on the ve- 
locity dispersion measured across the arms of the X in 
the V — 1 emission (Fig. [3]). We adopt the following 
relative (approximately solar) abundances: [SiO/H2] = 
10-4, PSiOI/P SiOl - 19.5, and psSiOj/poSiO] ~ 30.5 
(jPenziasI Il981f )). The E instein co efficients are calcu- 
lated from the data of iTipping fc Chackcrian (1981), 
while the collisi o nal ra te coefficients are taken from 
iBieniek fc Greml (fl983al lH). 

For simplicity, we assume that the velocity field is 
spherically-symmetric at any point in the medium, i.e. 
the logarithmic velocity gradient d{ln V)/dr — 1.0. 
The photon escape probability is then given hy P — 
(1 — e-'^=)/Ts, where the Sobolev optical depth is Ts — 
(c^/ Sit I'll) Aui[nu-nigu/gi]/{dV/dr), where dV/dr is the 
velocity gradient and and n; are populations per sub- 
level. Calculated population inversions are represented 
in terms of the gain coefficient 7 = An/Snts, where 
An is the population inversion and ts is the radiative 
lifetime of the transition. 

The physical conditions required for the excitation of 
different maser species are summarized in Table [5] The 
^^SiO v = OJ— 1 — transition is inverted at relatively 
low density (Nua < 10^ cm"'^), and it is never inverted in 
a strong, hot radiation field {W < 0.001), consistent with 
an inner radius of 200-1000 AU (Sect. [S]). It occurs at a 
wide range of kinetic temperatures but it is optimized at 
Tk < 1200 K. 

The ^^SiO u==l,2J=l — transitions are never 
inverted for Nh., < 10^ cm^'^ and they favor densi- 
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1) and 10^ 
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{v — 2). These results explain the spatial separation be- 
tween the ground (> 100 AU) and excited vibrational 
states (< 100 AU). The v — 2 transition is optimized 
at higher temperatures {Tk ^ 2000 K) and it is more 
strongly inverted in a strong, hot radiation field than 
the V — I transition, which is apparently quenched 
for W> 0.01. This is in agreement with the finding 
that V — 2 mascrs tend to lie closer to the star than 
V = I masers, although there is substantial spatial over- 
lap (Sect.[3l). 

The 29SiO and 3°SiO v = OJ=l-0 transitions are 
inverted across a broad range of parameter space. Their 
emission is optimized for conditions similar to ^^SiO v — 
1, 2 (Nh2 = 10« - 10" cm-3 and Tk = 1000 - 2000 K), 
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TABLE 4 

Input parameters in the LVG model. 





T,(2) 








x[SiO/H2](5) 


28SiO:29SiO:30SiO(5) 


dV/rfr(*') 


(K) 


(AU) 


(AU) 


(cm ^) 


(K) 






(km s-i AU-i) 


10* 


0.04-5 


10-1000 10-1 - 10-8 


104 - 1012 


600-2400 


lO-'i 


1:1/19.5:1/30.5 


0.5 



Note. — 

is the blackbody temperature of the radiative field 
'■^^ -Rrad is the radius of the radiative source. 
^"^^ -Rmas is the distance to the maser region, 
w is the radiation geometric dilution factor. 
x[SiO/H2] and ^^SiOi^^SiO 



SiO are the (solar) fractional abundances of SiO and isotopologucs. 

<«) dv/dr is the velocity gradient as measured by the velocity dispersion o^v^ <mW^c emission arises from compact Compo- 
nents with a spatial and velocity distribution gener- 
ally similar to the ^^SiO vibrationally excited emis- 
sion and traces material within 100 AU of Source I. 



TABLE 5 

Optimal LVG excitation conditions of 
SiO masers. 



Species 




W 












(cm ) 


(K) 


28SiO V = 





0.0001 


< 10^ 


< 1200 


28SiO V = 


1 


< 0.01 


108 - lOlO 


> 1500 


28SiO V = 


2 


0.01 


10» - 10" 


> 2000 


29SiO V = 





0.01 


108 - 10" 


> 1500 


30SiO V = 





0.01 


108 - 10" 


> 1500 



but can also occur more weakly under conditions con- 
ducive to ^^SiO V = inversion (Nhj < 10^ cm~^). 
These elements are in agreement with our finding that 
strong isotopic maser emission is excited within 100 AU 
from Source I as is ^*SiO v = 1,2 maser emission, while 
extended weak isotopic emission may occur at larger dis- 
tances from Source I in more extended structures, such 
as along the NE-SW outflow traced by the ^^SiO u = 
emission (Sect. |4T|) . 

In conclusion, our LVG radiative pumping model 
explains the bulk SiO maser emission characteristics 
of the five maser transitions observed toward Orion 
BN/KL. However, although the large-scale spatial and 
velocity distribution is similar for ■^^SiO v — 1,2 and 
^^SiO/^°SiO V — maser centroids, our data in- 
dicate differences in the small-scale spatial distribu- 
tions and a certain degree of anticorrelation in in- 
tensity. In order to take into account these effects, 
future work would need to incorporate non-local ra- 
diative transfer and possibility of line overlap among 
the isotopologuc s in the modeling of the exc i tation 
of SiO masers ( Gonzalez- Alfonso &: Cernicharol Il997l : 



iHerpin fc Baudry 



2000). 



5. SUMMARY 

We have used the VLA at 7 mm wavelength to image 
five rotational transitions from three SiO isotopologucs 
towards Orion BN/KL: ^^SiO w = 0, 1,2 and ^SSiO and 
30SiO ^; = (J = 1 - 0). We have imaged for the first 
time the ^^SiO and '^"SiO v ~ maser emission with an 
angular resolution ^0''2. 

The main findings of this paper are listed below. 

• ^*SiO V — maser emission traces a bipolar outflow 
driven by Source I extended over ~ 700 AU along 
a NE-SW axis. 



With resolution of 0;'2, ~ 70% of the isotopic 
emission appears to be resolved. The missing flux 
may arise in the outflow traced at larger radii by 
^*SiO V = emission. 



LVG radiative transfer calculations predict simi- 
lar temperatures and densities for the optimum 
maser excitation of both ^^SiO/'^°SiO v = and 
28SiO v = l,2 masers (T = 1000 - 2000 K, n^^ - 



10 



lOitl 



cm ), significantly higher than what is re- 



quired for ^«SiO w = (T ~ 600 - 1200 K, nn^ ~ 
10^^^ cm~^). However, weak emission from the 
isotopologucs cannot be excluded at the lower tem- 
peratures and densities that support ^'^SiO v = 0. 



• Although the emission distributions for different 
isotopologucs are generally similar (e.g., in radius), 
the peaks do not appear to overlap in detail (al- 
though blending effects may contribute to the ap- 
parent offsets owing to the coarse interferometer 
beam). 

The apparent positional offset of maser features from 
different transitions and/or isotopologucs suggests that 
non-local line overlaps should be included in the mod- 
eling of the excitation of SiO (isotopic) masers. Fol- 
low up modeling will use a non-local spherical geome- 
try radiative transfer code with line overlap. Higher- 
resolution and higher sensitivity simultaneous observa- 
tions of multiple maser lines/species with future facil- 
ities like EVLA and ALMA will be very important to 
make progress in the understanding of the gas dynamics 
around Source I and to better constrain physical condi- 
tions around Source I on small-scales. 



The data presented here were obtained under the VLA 
programs AG575 and AG776. This material is based 
upon work supported by the National Science Founda- 
tion under Grant No. NSF AST 0507478. 
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